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^ : Abstract 



We carry out a detailed study of exclusive radiative rare B-decays in the framework of 
the QCD sum rules on the light cone, which combines the traditional QCD sum rule 
Q^l technique with the description of final state vector mesons in terms of the light-cone 

1-Ch ' wave functions of increasing twist. Our calculation is restricted to the leading twist-two 

operators. The decays considered are: B^^ —* K* + 7, Bu^ ~^ P + 1, Bd ^ ou + •y 

k^ ', and the corresponding decays of the Bs mesons, 5^ — > + 7 and Bs -^ K* + 7. Based 

^ I on our estimate of the transition form factor Ff^^*{0) = 0.32 ± 0.05, we find for the 

branching ratio BR{B — > K* + 7) = (4.8 ± 1.5) x 10~^, which is in agreement with the 
observed value of (4.5±1.5±0.9) x 10~^ measured by the CLEO collaboration. We present 
detailed estimates for the ratios of the radiative decay form factors, which are then used 
to predict the rates for the exclusive radiative 5-decays listed above. This in principle 
allows the extraction of the CKM matrix element \Vtd\ from the penguin-dominated CKM- 
suppressed radiative decays when they are measured. We give a detailed discussion of 
the dependence of the radiative transition form factors on the 6-quark mass and on the 
momentum transfer, as well as their interrelation with the CKM-suppressed semileptonic 
decay form factors in B ^ p + i + ug, which we also calculate in our approach. 
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1 Introduction 

Theoretical interest in rare 5-decays lies in the first place in their potential role as preci- 
sion tests of the quark sector of the Standard Model (SM). Interpreted within this frame- 
work, the eventual measurements of these decays will provide quantitative information 
about the top quark mass and more importantly about the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements, Vtd, Vts and Vtb- In particular, the short- distance contribution to 
the CKM-suppressed rare -B-decays b —^ d + •y and 6 — > d + i^i~ directly measures Vtd- 
Together with improved measurements of the CKM matrix elements \Vcb\ and \Vub\, this 
will determine the CKM unitarity triangle resulting from the constraint J2i VidVil = 0, 
and thus, pin down the CP-violating phases in the Standard Model. These non-trivial 
constraints on the CKM unitarity and the importance of rare B-decays in this respect 
cannot be overemphasized. At the same time, rare S-decays induced by flavour changing 
neutral currents (FCNC) have the potential of providing one of the early hints for non-SM 
physics. It is, therefore, imperative to get estimates of these FCNC decays in the SM as 
reliable as possible and carry out an experimental physics programme sensitive to rare 
B-decays. 

The experimental searches for the FCNC -B-decays has already provided first divi- 
dends. The recent CLEO observation [^ of the rare decay mode B -^ K* + 7 having a 
combined branching ratio BR{B —>■ K* + 7) = (4.5 ± 1.5 ± 0.9) x 10~^ and an improved 
upper limit on the inclusive branching ratio BR{B ^ X^ + 7) < 5.4 x 10~^ (95% C.L.) 
0] have been analysed in the SM context, as well as in several SM extensions. Within the 
SM, the resulting experimental measurements have been interpreted in terms of upper 
and lower limits on the CKM matrix element ratio 0.50 < |l^s|/|Vcfe| < 1.67 (at 95% C.L.) 

i- 

On general grounds, the most quantitative tests of FCNC processes in B-decays are 

the ones involving inclusive measurements such as B — >■ X^^s + 7 and the corresponding 

semileptonic and purely leptonic processes such as B — >• X^ ,, + i^£~ and B^^ ,, — ^ £+£". 

Thus, the ratio of inclusive decay rates such as BR{B -^ Xd + 'j)/BR{B -^ Xs + 7) would 

provide one of the cleanest determinations of |Vjrf|/|\4s| [^. However, inclusice branching 

ratios are rather difficult to measure. Exclusive decays such asB^p + 7, B^co' + 7 and 

the semileptonic decays B -^ (vr, p, u)i^£~ are probably easier to measure, given a large 

number of B hadrons and a good electromagnetic detector and particle identification. The 

interpretation of these and the already measured decay B — > K* + 7 necessarily needs 

reliable estimates of the decay form factors. The decay B -^ K* + 7 has received quite a 

bit of theoretical attention p|-p!5|. 

The aim of this paper is to estimate a number of exclusive radiative decays. These 

include the CKM-allowed decays, B„ ^ -^ -ft'*+7 and B^ -^ 0+7, and the CKM-suppressed 

decays, B^^ ^ p + j, B^ ^ uj + j and Bg — > K* + 7. While the effective Hamiltonian 

approach allows to include short- distance QCD corrections at scales /i^ > ml, estimates 

of the hadronic matrix elements of the relevant operators necessarily require some non- 

perturbative technique. We use here QCD sum rules, which have proved to be a powerful 

tool for such purposes, and in their classical form have been introduced in the pioneering 



papers by Shifman, Vainshtein and Zakharov |16]. In this paper, we use a modification 



of the QCD sum rule approach, which we refer to as QCD sum rules on the light cone. 



This technique has been developed originally for light quark systems in [T7, |TB[. The 
application for the heavy meson decays was first suggested and studied in |jl9|. 

In this approach, the ideas of duality and matching between the parton and hadron 
descriptions, intrinsic to the QCD sum rule framework, are combined with the specific 



techniques used in the studies of the hard exclusive processes in QCD [^ ^. In contrast 



to the standard sum rule method for the evaluation of form factors [£^) |2^ ? iii which the 
hadrons in the initial and final state are treated in a symmetric way, we resort to the 
QCD sum rule treatment of the initial 5-meson only, and describe the outgoing light 
vector meson by the set of its wave functions of increasing twist. Hence the operator 
product expansion in the light-cone approach is governed by the twist of the operators 
rather than by their dimension, and the vacuum expectation values of local operators are 
replaced by the light-cone hadron wave functions. A physical motivation of the light-cone 
sum rule approach in the present context is an obvious asymmetry of the participating 
heavy initial and light final state mesons. The advantage of this formulation is that it 
allows to incorporate additional information about high-energy asymptotics of correlation 
functions in QCD, which is accumulated in the wave functions. The high energy behaviour 
of these functions is related to the (approximate) conformal invariance of QCD, and many 
properties and results, following from this approximate invariance and obtained previously, 
can be advantageously used in the present context as well. Existing applications of the 
light-cone sum rules include the calculations of the amplitudes for the radiative decay 



S — > p7 [0, nucleon magnetic moments []TB|, the strong couplings gj^NN and gp^^T, |T^ 



the semileptonic S-meson |jT9|, ^ and D-meson decay amplitudes |2^. In all these cases 
the results have been encouraging. 

The principal theoretical result of this paper is a sum rule for the electromagnetic 
penguin form factor appearing in the decays B ^ V + 'y, where V^ is a vector meson. We 
derive the corresponding sum rule in the traditional QCD sum rule approach also and 
comment on similar existing sum rules in the literature 0-|@]. We argue that the sum 
rules derived in the light-cone approach are more reliable and we calculate the radiative 
transition form factors for a number of exclusive radiative 5-decays in this framework. 
The main numerical results of our work can be summarized as follows: 

Ff-^^' = 0.32 ±0.05 pB-^ip^-^) = 0.24 ±0.04 , . 

pB^^d, ^ 0.29 ±0.05 pBs^K' ^ 0.20 ±0.04 ^ ' 

where the form factor Fi is defined below (eq. (^)). The above estimate for the decay 
form factor Ff^^ can be combined with the QCD-improved inclusive radiative branching 
ratio BR{B ^ X, ± 7) = (3.0 ± 1.2) x 10"^ § to yield BR{B ^ K* + -f) = (4.8 ± 1.5) x 
10~^. This agrees fairly well with the corresponding branching ratio measured by the 
CLEO collaboration, posted as (4.5 ± 1.5 ± 0.9) x 10"^ [|T|. Together with the eventual 
measurements of the CKM-suppressed radiative decays B^^ — >p±7, B^ ^ ou + •y and 
Bg — > K* ± 7, the remaining form factors in (||) can then be used to determine the CKM 
parameters. 

While the main thrust of the paper is on providing theoretical estimates of the ex- 
clusive radiative S-decay rates, we also discuss a number of other, related, theoretical 
and phenomenological issues. In particular, we discuss in detail the heavy-quark-mass 
dependence of the radiative 5-decay form factors, and derive the sum rules in the heavy 



quark limit. A light-cone QCD sum rule for the semileptonic decay B ^ p + iui is also 
derived and the relation between radiative and semileptonic form factors is studied. 

The paper is organized as follows. In section 2, we review the effective Hamiltonian 
for radiative 5-decays, dominated by the electromagnetic penguins, and we work out 
explicitly the dependence of their rates on the CKM matrix element parameters. A 
number of relations involving exclusive and inclusive radiative i?-decays is given in this 
section. In section 3, we describe the method of the QCD sum rules on the light cone and 
derive our general sum rule for the form factors of the electromagnetic penguin operators 
entering in the decays B ^ V + ^. The wave functions of the vector mesons which are 
needed to evaluate the sum rule are discussed in section 4, and the numerical estimates 
of the form factors are presented in section 5. Section 6 includes the discussion of the 
6-quark mass dependence of the form factors, and the heavy quark limit. In section 7 
we discuss the g^ dependence and the interrelation of radiative and semileptonic -B-decay 
form factors in the QCD sum rule approach. Section 8 contains a summary and some 
concluding remarks, and the impact of measuring the decays 5^ — > u; + 7, Bu^ — ;> p + 7 
and Bg -^ K* + 7 in determining the CKM matrix element \Vtd\ is underlined. The 
traditional three-point QCD sum rule for the radiative -B-decay B —>■ K* + 7 is derived 
and discussed in the appendix. 

2 Effective Hamiltonian for inclusive and 
exclusive radiative B-decays 

In this section we review the theoretical framework within the Standard Model of elec- 
troweak interactions for inclusive and exclusive radiative 5-decays dominated by the elec- 
tromagnetic penguins. The inclusive decays are grouped as -B — > Xf + 7, where we use 
the flavour of the light quark / = s, c? in the transition b ^ f to characterize the hadronic 
system recoiling against the photon. Including lowest-order QCD corrections and gluon 
bremsstrahlung, these decays are described at the parton level by the transitions 6—^/7 
and b -^ f'^g- In calculating the inclusive decay widths, we shall follow the work reported 
in [^, ^, |TT|. The exclusive decays that are the principal concern of this work are grouped 
in an analogous way into b -^ s transitions: 

. 5„^ir* + 7, 5rf^ir* + 7, 

• 5, ^ + 7, 

which we shall also term CKM-allowed, according to the dominant CKM matrix element 
dependence of their decay rates, and b ^ d transitions: 

• 5d ^ p + 7, Bd^uJ + i, Bu^ p + 'j, 

• B,^K* + ^, 

which will be called CKM-suppressed. 

The framework to incorporate (perturbative) short- distance QCD corrections in a 
systematic way is that of an effective low energy theory with five quarks. It is obtained 



by integrating out the heavier degrees of freedom, i.e. the top quark and W^ bosons. 
Since the running of a^ between mj (present estimates rrit = 164 ± 27 GeV [^) and mw 



is not very significant, and since mw/rnt would not be a good expansion parameter, it is 
a reasonable approximation to integrate out both the top quark and the W^ at the same 
scale. 

Before using the unitarity properties of the CKM matrix, the effective Hamiltonian 
relevant for the processes 6—^/7 and b —>■ f'jg has the form 

AC ( ^ 

^S7'^ = -^\v,,v:fY:cAii)o,{i^) 

+ K.K}Ec;(/i)o;(/.) + K.K*/E^»o» , (2) 

where Gp is the Fermi coupling constant and j runs through a complete set of operators 
with dimension up to six; the Cj(/i) are their Wilson coefficients evaluated at the scale /i. 
We recall here that O^ and O2 represent the colour-singlet and colour-octet four- 
fermion operators, respectively, obtained from the SM charged current Lagrangian written 
in the charge retention form; in particular, 

O2 = {cLa1^bLa){fL(BlfiCLt^) and O'^ = {uLal^bLa){fLl3lfiULl3) , (3) 

where a and (3 are SU(3) colour indices. The remaining operators are the same for all 
three CKM prefactors, i.e. O' = O" = Oj for j ^ 1,2. At tree level, the only contribution 
to b —>■ f'-f comes from the magnetic moment operator 

Or = yI^ f^y^" {m,R + rrifL) b F,, , (4) 

where L,R = (l=P75)/2, and e is the QED coupling constant. The four-fermion operators 
O3, ...,0e and the QCD magnetic moment operator Og, which is the gluonic counterpart 
of O7, arise from penguin diagrams in the full theory (before integrating out W and t). 
These operators enter indirectly in 5 — >■ /7 decays due to operator mixing and through 
(virtual and bremsstrahlung) gluon corrections. 

Taking into account the unitarity of the CKM matrix, only two terms in the combi- 
nations Aj = VibV*. — or in the combinations ^j = VibV*^ in the case of CKM-suppressed 
transitions — are independent. For b —>■ s transitions, one finds |Au| <^ |Ac|, |At|; therefore, 
when neglecting terms proportional to A^, the effective Hamiltonian (§) can be brought 
into a form proportional to A^ ~ —Ac: 



<r^=-^A,EQ(/^)o,(/x) . (5) 



For b ^ d transitions, where C,u, ^c and C,t are all of the same order of magnitude, it 
is most convenient to choose C,u and C,c as independent CKM factors during the matching 
for the Wilson coefficients at /x = mw and during the renormalization group evolution to 
fi ~ nib- Since terms of order 0{ml^/m^r) and 0{ml/mly) (or O^ml^/mf) and 0{ml/mf)) 



are thereby neglected, the terms proportional to ^„ and ^c are multiplied by just the same 
coefficient functions, i.e. Cj(/x) = C'-{^) (although the corresponding operators are, of 
course, different, see (^). Finally, exploiting again ^u + ic = —C,t, one ends up with 

where the Wilson coefficients Cj{fi) are precisely the same functions as in (|). Details of 
the full operator basis, the matching of the Wilson coefficients Cj at /i ~ mw, and the 
complete leading-logarithmic renormalization group evolution to /x <^ niw can be found 
in the literature 



In estimates of the inclusive branching ratio for BR{B — i> X^ + 7) one has to include 
the contribution of the QCD bremsstrahlung process 6 — *> s + 7 + (7 and the virtual 



corrections to 6 — *> s + 7, both calculated in 0{aas) in ref. [|ri|]. Expressed in terms of 



the inclusive semileptonic branching ratio BR{B -^ Xluf), one finds: 
BR(B -. X. + ,) = ,-—-^-^-^^^^-^^^^j^-^--^ X BB.(B -. XU„) . (7) 

where Xt = ml/m^ indicates the explicit mf-dependence in Ci{xt,mf,) = Ci^nif,), and 
g{r) = 1 — 8r^ + 8r^ _ 7-8 _ 24r^ ln(r) is the phase-space function for r(6 ^ c + iui). The 
function /(r) accounts for QCD corrections to the semileptonic decay and can be found. 



for example, in ref. |32]. It is a slowly varying function of r, and for a typical quark mass 
ratio of r = 0.35 ± 0.05, it has the value /(r) = 2.37 =F 0.13. The contributions from 
the decays b —>■ u + d, vp have been neglected in the denominator in (|^) since they are 
numerically inessential (|V^b| <^ iKftD- For the semileptonic branching ratio, the measured 
value BR{B -^ Xlue) ~ 11% will be used. 

The inclusive decay width for B ^ Xs + '~f is dominantly contributed by the magnetic 
moment term Cj{xt-, ^i)0-j{^) , hence the rationale of factoring out this coefficient in the 
expression for BR{B — > X^ + 7) in Eq. (|^) above. Including 0{as) corrections brings 
to the fore other operators with their specific Wilson coefficients. The effect of these 
additional terms can be expressed in terms of the function K{xt, rub), which is a i^-factor 
in the sense of QCD corrections. The function K{xt, rub) has been computed in |Q taking 
into account the dominant corrections from C2 and Cg (the coefficients of other operators 
are considerably smaller [^). The resulting branching ratio is0 0: 

Si?(5^X, + 7) = (3.0±0.5) X 10"^ (8) 

for rrit in the range 100 GeV < rrit < 200 GeV. This is to be contrasted with the present 
upper limit on the exclusive radiative decay, BR{B — *■ Xs-l-7) < 5.4 x 10~^ (at 90% C.L.) 

i- 

It has been argued in ref. that the scale dependence of Cji^Xt, /i) is very pronounced 
in the presently available leading-logarithmic approximation. This inherent uncertainty 
of the present theoretical framework has to be borne in mind in a quantitative discussion 

■f obtained by setting the CKM matrix element ratio |Atp/|Vc6p to 1, according to our present knowl- 
edge of the CKM matrix (readily seen, e.g., in the Wolfenstein representation ||27[|). 



of the QCD-improved decay rates. In ref. this uncertainty has been estimated by 
varying the scale parameter n in the range mf)/2 < fi < 2mh and introduces an additional 
theoretical error of ±1 x 10~^ in the above estimates for BR{B —>■ Xg + j). 

The case of inclusive B -^ Xa + 7 decays is somewhat more complicated with respect 
to the factorization of the CKM parameters. When evaluating the one-loop matrix ele- 
ments for the parton process b —>■ s'jg, at least the four-quark operators O2 and O2 should 
be kept since their coefficients are of order unity (they contribute in penguin diagrams 
emitting both a gluon and a photon). These contributions involve a non-trivial depen- 
dence on the CKM matrix elements through terms proportional to ^„ and ^c, see eq. (^. 
Explicit expressions for this dependence on the p and r) parameters in the Wolfenstein 
representation of the CKM matrix have been evaluated in ref. [|]. 

In the following we shall focus our attention on exclusive decays, like B^^ -^ K* + 7 
and -Bs — >■ + 7, and the corresponding CKM-suppressed modes listed at the beginning 
of this section. At tree level only the magnetic moment operator Of of eq. @ contributes 
to the transition amplitudes. For a generic radiative decay B ^ V + "^ one defines a 
transition form factor Fi{q^) as: 

{V, \\fa^,q^h\B) = ze^,,.e*W^p^p^2Ff ^^(g2) , (9) 

where ^ is a vector meson [V = p,uj,K* or (p) with the polarization vector e^'^-'; and 
B is the generic 5-meson B^.B^ or Bg. The vectors p^, Vv and q = Pb — Pv denote 
the four-momenta of the initial 5-meson and the outgoing vector meson and photon, 
respectively. 

In (H) it is understood that the operator is evaluated at the scale p = rrif,, and all large 
logarithms, \ia{mw/fnb) and ln{mt/mb), are included in the coefficient function C7(/x = 
mb). The 6-quark mass, which has been factored out, should be identified with the pole 
mass, although the complete two-loop treatment of the coefficient function is needed to 
make this identification meaningful. 

In evaluating the hadronic matrix elements, one may consider the 6-quark mass as 
a large parameter, and try to collect logarithms, corresponding to the so-called hybrid 
anomalous dimension |R3[. At zero-recoil, g^^^^ = (m^ — rnvY — 19 GeV^, this treatment 



is simple and the answer is obtained in the framework of the heavy quark effective theory. 
However, in this case the extrapolation to the physical point q^ = introduces a large 
uncertainty. We take a different approach, and use the QCD sum rules to calculate the 
form factors directly at the physical point q^ = and for the finite value of the 6-quark 
mass. The price to pay is that the treatment of hybrid logarithms becomes complicated, 
and we shall ignore them in this paper. 

With the above definition, the exclusive decay widths are given by {B = Bu or Bd)'- 

nB - K^ + 7) = ;^G||A,| Vf-^'*(0)pC,(x„ mbf{ml + ^2) ("^|-<0^ ^ ^^q) 

and the analogous expression for T{Bs -^ + 7). The branching ratios of these exclusive 
decays can again be written in terms of the inclusive semileptonic branching ratio: 

Hn(H K^ I -^,1 -fi" l^^l' |C7(x„m,)nFf-^*(0)P (l-m^.Vm|f 

^ ^^^ 7:\V,b? g{mJmb){l-2/?,^f{mJmb)) [l-ml/mlf ^ '°^- 

(11) 



A good quantity to test the model dependence of the form factors for the exclusive 
decay is the ratio of the exclusive-to-inclusive radiative decay widths {B = B^ or Bd): 

..^*/^ ^ = V{B^K* + ^) ^ (l-m,,,Vm|)3 m| |Ff'-^'(0)p 

Note that K{xt, rrif,) ~ 0.83 is almost independent of rrit (for the range 100 GeV < 
rrit < 200 GeV) 0. The exclusive-to- inclusive ratio involving S^-decays is defined in an 
analogous way. 

Since the same short-distance-corrected coefficient function Cj{mij) enters in the Hamil- 
tonian for CKM-allowed (|^) and CKM-suppressed (|^) modes, the QCD scaling is identical 
for the two-body decays 6 ^ s + 7 and 6 ^ (i + 7, and it does not affect the ratio of the 
decay widths V{h ^> d + '~^) /V{h -^ ■5 + 7)- The same applies for the exclusive decays such 
as i? — > p + 7 and B -^ K* + 7, in which case the CKM factors factorize in the decay 
amplitudes. 

From this observation a number of relations between the exclusive decay rates follow 
in the Standard Model [Q. This is exemplified by the decay rates for B^^d ~^ P + 7 and 

Bu,d^K* + r. 



r(5„,rf^ir*+7) \Vts\'\F^-''' 
where (^u,d is a phase-space factor: 



;^u,d. (13) 



(ml + ml) {nil^^,-mlf 



The ratio (|13D depends only on the CKM matrix elements and the ratio of form factors, 
while it is independent of the top quark mass (and of the renormalization scale fi). 

Note that the decay width T{Bd — > u; + 7) is expected to be equal to the decay width 
T{Bd — >■ p + 7), apart from the minor difference in the phase-space factors $. This follows 
from the assumption that the quark wave functions for uj and p are described by the 
isoscalar and isovector combinations, \uj) = l/^/2{uu + dd) and \p) = l/\/2{uu — dd), 
respectively. 

The CKM-suppressed exclusive decay width T{Bs — > K* + 7) can be related to the 
CKM-allowed decay width T[Bd -^ K* -|- 7) by a relation similar to the one given in 
eq. (p!3|) . However, one expects a substantial difference in the form factors F^ '^^ (0) and 
Fi '"^ (0) due to the exchange of the roles of s and d quarks in the wave function of the 
K* in the two decay modes. 

3 QCD sum rules on the light cone 

The aim of this and the following two sections is to calculate transition form factors, 
governing the radiative -B-decays B ^ V + ■y, as defined in (^. Our approach is very 
close to the calculation of the semileptonic B -^ -kcv form factor in |l^, ^. 
To derive the sum rule, we consider the correlation function 

I fdxe'''^V{p,X)\T{^p{x)a^,q''b{x)mzj,ij{0)}\0)=ze^,p,e<^>q'p''T{{p + qf) (15) 



at q^ = 0, p^ = rriy^ and at Euclidean m1 — {p+qY of order several GeV^. Hereafter we use 
^ as a generic notation for the field of the light quark. Writing down the dispersion relation 
in {p + qYi we can separate the contribution of the 5-nieson as the pole contribution to 
the invariant function T((p + qY): 

s9s /b^b 2Fi(0) , , 

T{{p + qf)= ^\^ — r^ + --- ' ^6 

mi, + niq m% - [p + qY 

where the dots stand for contributions of higher-mass resonances and the continuum. The 
5-meson decay constant is defined in the usual way, 

(0|V^7m75&|5(p))=^Pm/b, (17) 

and mBi frif, and rriq are the 5-meson, 6-quark and light-quark masses, respectively. 

The virtuality of the heavy quark in the correlation function under consideration is 
large, of order ml — {p + qY, and one can use the perturbative expansion of its propagator 
in the external field of slowly varying fluctuations inside the vector meson. The leading 
contribution corresponds to the diagram shown in Fig. la, and equals 



d^^""" / — — e-^"" ,^" , /\/(p,A)|T{^(a;)a^.K+ ^)275^(0)}|0). (18) 



^igx I dk jfc3, qu 

To the leading-twist accuracy, taking into account gluon corrections like the one in Fig. lb 
produces the path-ordered gauge factor in between the remaining quark fields in (pTSD, 
which we do not show for brevity. 

Thus, in general we are left with matrix elements of gauge-invariant non-local oper- 
ators, sandwiched in between the vacuum and the meson state. These matrix elements 
define the light-cone meson wave functions, which have received a lot of attention in the 
past decade. Following Chernyak and Zhitnitsky |^, we define 



{o\m^,uH^)v{p,X)) 


= ^(ei'^P 


u-e^^^P,)f^l' 


due 


'■"P^0^(u 


/.^). 


(19) 


Likewise, 
















{0m0h,^{x)V{p,X)) = 


(e^x) 


fvmv / du 
Jo 


e-»^>ll(M,/i2) 








\ 


(eWx)\ 
^^ ipx) ) 


fvmv 

Jo 


t/we-^^P" 


9i\u, 


/i'),(20) 


(0^(0)7^75^(x)F(p,A)) = 


1 


Jo 


due 


-^-^^g^_^\ 


u,^?). 


(21) 



The functions 0_l(u, /x^) and 0||(u, /i^) give the leading- twist distributions in the fraction of 

total momentum carried by the quark in transversely and longitudinally polarized mesons, 

respectively. The functions g^ {u,^"^) and g]^ {u,^'^) are discussed in detail below. The 

normalization is chosen in such a way that for all four distributions / = 0_l, 0||,f7_L ,(?_]" , 

we have: 

/ du f{u) = 1. 
Jo 



In the matrix elements of non-local operators on the l.h.s. of (|T^-(^) the separations 
are assumed to be light-like, i.e. x^ = 0. Regularization of UV divergences that arise 
in the process of the extraction of the leading x^ — > behaviour produces a non-trivial 
scale dependence of the wave functions, which can be found by renormalization group 
methods [^, ^. The scale in (|18]) is fixed by the actual light cone separation, ij? ~ 



X 



~ m'. 



-{p + qf. 



Putting eqs. (p!^)-(|2T|) together, we obtain 



T{{p + qy 



du 



ml + uuiTLy — u{p + qy 



mbfv<P±{u) + umvfv9±\u) 



+ -^myfyg^r^ 



+ - du 



mt 



unit 



4 Jo {ml + uurriy — u{j) + qY) 



;mvfv9± (u). 



H, 



:22i 



where u = 1 — u. The expression in (p^ has the form of a dispersion integral in [p + g)^, 
which can be made explicit by introducing the squared mass of the intermediate state 
s = m^l/u + um^Y as the integration variable, instead of the Feynman parameter u. The 
basic assumption of the QCD sum rule approach is that the contribution of the i?-meson 
corresponds in this dispersion integral to the contribution of intermediate states with 
masses smaller than a certain threshold s < so (duality interval). Making the Borel 
transformation l/(s — (p + q)"^) — > exp(— s/t) and equating the result to the 5-meson 
contribution in fllGf), we arrive at the sum rule 
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mvfv9±{u,iJ, =t) 



(23) 



which should be satisfied for values of the Borel parameter t of order several GeV^. To 
the leading logarithmic accuracy, the scale in the wave functions coincides with the Borel 
parameter. Principal input in this sum rule are the vector meson wave functions, which 
contain non-trivial information about the dynamics at large distances, and which we are 
going to discuss now. 

4 Wave functions of the vector mesons 



It is known that the decomposition of the leading-twist meson wave functions in terms 
of conformal invariant operators allows one to diagonalize the mixing matrix at one-loop 
order. Thus, the (approximate) conformal invariance of QCD implies that the coefficients 
in the expansion of leading-twist wave functions in the series of Gegenbauer polynomials 



37| are renormalized multiplicatively to the leading logarithmic accuracy. In our case 
0j_(M,/i) = 6m(1 — m) 



1 + ai(/i)e + a2(/i) ( ^' - 5) + «3(/i) (3^' 



e + 



an(/^) = a„(/io)f^^j , (24) 

where we have introduced the shorthand notation ^ = 2m— 1. Here b = {ll/3)Nc—{2/3)nf. 



The anomalous dimensions turn out to be 138 



7n = C^fn-4£l/A (25) 



where Cp = {N^ — l)/{2Nc). The coefficients in the Gegenbauer expansion ( p4|) at a 
low scale, a„(yUo), should be determined by a certain non-perturbative approach, or taken 
from experiment. At present, most of the existing information comes from the QCD sum 
rules. Following |^, we use the model wave functions for p, K* and mesons at the 



scale /^o = 1 GeV , corresponding to the following choice of the parameters in (|2 

4"^ = -1.25, 
oP = 0.75, aP = -2, aP = 0.75, 

a^f = -2.5. (26) 

Note that in the case of p and mesons only 02 is non-zero, and the wave function is 
symmetric under the interchange u ^^ 1—u. In the case of a ii'*-meson, it is necessary to 
specify that the variable u is the fraction of the total momentum carried by the strange 
quark. The decay constants appearing in (|I^) are [^: 

ff = 200 MeV, 
/^. = 210 MeV, 
/^ = 230 MeV. (27) 

As mentioned above, the correct normalization point in the sum rules is of the order of 
the typical Borel parameter, which for 5-meson decays is of order /i^ ~ m^ — ml ~ 5 
GeV^. Using Ajjg = 225 MeV in the renormalization group rescaling factors given in 
(1^), we obtain at this scale 

4^) = -0.85, 
af*^ = 0.57, 4^^*^ = -1-35, 4^'^ = 0-46, 

ai^^ = -1.7. (28) 

The wave functions corresponding to the above parametrization are shown in Fig. 2. As 
a general effect of the rescaling, the wave functions become somewhat wider and closer 
to the asymptotic expression 6m(1 — u). In the case of the ii'*-meson, the wave function 
becomes both wider and more symmetric, so that scaling violation affects the region 
M < 1/2 only ( see Fig. 2b). The integration over u in the sum rule in ( ^31) is restricted 
to the interval of rather large momentum fractions, carried by the light quark involved 
in the electromagnetic penguin operator. For realistic values of parameters, see below. 



one has u > 0.65-0.7. According to the analysis in [RH], the K* wave function turns out 



to be in this region very close to the asymptotic expression 6^(1 — u). In the case of 
the decay Bs —>■ K*'~f, the role of the strange and non-strange quarks in the i^*-nieson 
is reversed, corresponding to the formal substitution u — *■ 1 — u in the wave function. 
Thus, the region of small values u < 0.3 — 0.35 in Fig. 2b becomes appropriate, where 
deviations from the asymptotic form are large and the effect of rescaling is important. 
This results in a significant difference between the form factors in the decays B^i -^ K* + '-f 
and Bs-* K* + 7. 

The conformal expansion of the wave functions g± (u, jj?) and g± (u, /i^) is somewhat 
more involved and can be obtained using the approach of ||39| . To this end we define new 
wave functions, which contain quarks with fixed spin projections on the light cone: 









(29) 



The wave functions g^^{u) and g^\u) can be expanded in terms of irreducible representa- 
tions of the collinear conformal group as pO|, ffTl p9| : 



g^^u) 

g^Hu) 



Ji 0(1,0) 



Ti 0(1,0), 



2(1 -n) i+c['pr'io+c'2'pr'io 



2u 



JT 0(0,1) 



iTD(0,l)/ 



i+c-^pr'io+c-^pr'io + 



(30) 



where PJ, '"^ (0 ^ire Jacobi polynomials [0. The normalization in (pOf ) follows from the 
standard definition of vector decay constants 



mi,i^\Vip,\)) = e'myfv. 



(31) 



The numerical values determined from QCD sum rules are |16, 36 



fp 
Ik* 



200 MeV , 
210 MeV , 
230 MeV . 



(32) 



Note that these couplings, to the accuracy of the existing QCD sum rule calculations, 
coincide with the couplings in (|27|). 

Neglecting S'f/(3)-breaking effects related to the difference of the quark and antiquark 
masses in the mesons, one has 



Cfc 



\fc JT 



-W^ = cu. 



(33) 



Finally, the coefficient Ci is actually fixed by the equations of motion, which allow one to 
reduce the matrix element 



^^C^i/po-C p 



(0|V'7m75(^/^.)V^|V^(p,A)) 
to the matrix element in (|3lD. Thus, we find A = —{l/2)fvmv ||4^ and 

ci = -l/2. 



(34) 



(35) 



On the other hand, one has the obvious relations: 



9^l\u) = -[gn(u) + g^Hu)] 



du 



g'l>{u) = 2[g^^{u) - g^^iu) 



(36) 



Using the identities [^ 

(1 + OPlf''\0 + (1 - OPt'''\0 = '^cl'\o 
(1 + 0Pi"''\0 - (1 - OPi^^'\i) = 2c}iUo 



d_ 



j2\^3/2 



1/2 



(1 - e)cr{o = -ik +i)ik+ 2)c;,i\io 



(37) 



where C^(0 ^^^ Gegenbauer polynomials, we arrive at the expansions: 



-.i^)i 



k=2n 



91' W = ^ (Cfc - Cfc_i)C^/ (0 , 

Cfc -Ck 



g^^{u) = 2(l-aE77:^^T^^^f(0 



(3^ 



The coefficients in front of each Gegenbauer polynomial come from the operators with 
two neighbouring conformal spins. According to a general result 0], the operators with 
different conformal spin do not mix under renormalization to leading logarithmic accuracy. 
This is a major simplification but does not guarantee multiplicative renormalization in the 
present case because of the existence of three-particle antiquark-quark-gluon operators 
with the same twist and conformal spin (see [pQp. 



Further insight in the structure of the wave functions g\_ {u 



(v), 



and g± (u) 



can be ob- 



tained by using the equations of motion. Note that both these wave functions correspond 
to transverse spin distributions, and involve precisely the same operators (apart from a 
missing 75 in the case of g]^ (u)) as the ones involved in the operator product expansion 
for the structure function g2{x,Q'^) of polarized deep inelastic lepton-nucleon scattering. 
The difference is indeed in the matrix elements - which involve the (polarized) nucleons 
with equal momenta in the latter case and the vacuum and the vector meson state in 
the present situation. Thus, similar to the case of the structure function g2, the wave 
functions g]^ (u) and g^^ (u) contain contributions coming from both operators of twist 
2 and twist 3. In close analogy to deep inelastic scattering, the twist-2 contributions to 
the "transverse" wave functions g_l , g^ can be expressed in terms of the leading- twist 
"longitudinal" wave function (j)\\{u, ij?) defined in 



(d), twist— 2/ \ 

gY' [u) 

d (a),twist-2/ X 



"" dv^-^ + t dv^-^ 



Jo V Ju V 



(39) 



which is the analogue of the Wandzura-Wilczek relations P3| between the structure func- 
tions g2 and gi. The remaining twist-3 contributions to g^ , g\^ can be written in terms 



of three-particle antiquark-gluon-quark wave functions of transversely polarized vector 
mesons, considered in |]4^ , p6| . The derivation of ( ^Of ) and the relations for twist-3 con- 
tributions will be given elsewhere. 

In this paper, we do not take into account contributions of twist 3 which come from 
the gluon-exchange diagram in Fig. lb, and to this accuracy we need to keep the twist-2 
contributions to the wave functions g]^ , g]j only. In anology with a similar situation in 
the case of the structure function g2, one should expect that corrections to the asymptotic 
wave function 0||(m) = 6m(1— m), coming from twist-2 operators involving gluons p6| are of 



the same order of magnitude as neglected contributions of twist 3. Thus, for consistency, 
we must neglect the gluon corrections to gl\ gl altogether, which amounts to putting 
to zero all the coefficients in (|30D except the first two, Cq = 1 and Ci = —1/2. Thus, we 
obtain, finally 

g^l\u) = ^(l-e)=6u{l-u), (40) 

which we use in the numerical analysis. It is easy to check that these expressions are 
exactly the ones that follow from (^) with the asymptotic expression for the longitudinal 
distribution amplitude 0|j(u) = Qu{l — u). The same expressions for the transverse wave 
functions in ( ^OD have been found in |^2|] by a different method. The identification of 
these contributions as being of twist 3 in ref. [H2| is, however, an error. 



5 Evaluation of the sum rules 

To complete the list of entries which appear in the sum rule (p3D we must specify the value 
of the pole 6-quark mass mi, and the continuum threshold Sq. We vary these parameters 
in the limits 

rrib = 4.6-4.8 GeV , 
sj,^) = 33-35 GeV^ , (41) 

and for the continuum threshold in the i^^-channel, we take 

sBs _ ,iJ ^ ^1^ _ ^1 ^ 1 GeV^, (42) 



making use of the present world average ms^ = 5373.2 ± 4.2 MeV [Q, ^, |4g . 

The value of the S-meson decay constant Jb has received a lot of attention recently. 
Within the QCD sum rule approach, its value is strongly correlated with the values of 
the 6-quark mass and the continuum threshold (see [^ and references cited therein). For 
consistency, we substitute the value of fs in (|23| ) by the square root of the corresponding 
QCD sum rule: 

JbJK -iml-ml)/t 
2 

mt 



^7r2 



flsdse (' <'^/\l - ml/s)^ - mb{qq)f,2^t - mb{qagGq)^2 

' (43) 




in which we discarded the radiative 0{as) corrections, as they are not taken into account 
in the sum rule (|23|) either, and we also neglected numerically insignificant contributions of 
the four-quark operators and of the gluon condensate.| The sum rule in (^) is evaluated 
at precisely the same values of nib and sq as in (pSf), and using the standard values of the 
vacuum condensates {qq) = —(240 MeV)^ and {qcrgGq) / {qq) = 0.8 GeV^ (at the low scale 
/Xq = 1 GeV^), which have been used in the QCD sum rule analysis of the wave functions 



in [^. The renormalization group evolution is given by 



2/(36) 



,"s(yUo), 
,«s(/^o), 



(qagGq),. = ^^ {qcrgGq)^2. (44) 



For the decay constant of the i?s-meson, we use the sum rule in (|^), and the ratio 

fBjfB, = 1-15 ±0.05, (45) 



which is in agreement with most existing estimates [g, ^ . 

The stability plots for the form factors as functions of the Borel parameter are shown 
in Fig. 3. The stability is in all cases good, and the variations of the continuum threshold 
within the limits (^) induce uncertainties in the values of form factors within less than 
3%. Using the region in the Borel parameter 6 GeV^ < t < 9 GeV^, we extract the values 
of form factors 

Ff^^' = 0.32 ±0.05 

pB^(p,i.) ^ 0.24 ±0.04 

pBs^^ = 0.29 ±0.05 

pB,~.K* ^ 0.20 ±0.04 (46) 

The given error mainly comes from the uncertainty in the 6-quark mass and the depen- 
dence on the Borel parameter. In addition, we have included an uncertainty of 20% for 
the parameters ( PUj) of the wave functions. This leads for all decays to an error of order 
±(0.01-0.015) for Fi. For all decays except Bg -^ K*'~f, the sum rule is dominated by the 
contribution of the wave function (f)j_, and the average value of the momentum fraction u 
under the integral in (p3[) is {u) ~ 0.8. 

The uncertainties in the values of input parameters tend to be reduced in the form 
factor ratios, see Fig. 4, and we obtain 

pB^K* =0-76 ±0.06. (47) 

§ The sum rule in (|4^) yields /b ^ 120-160 MeV, which is lower than the value preferred at present, 
typically fs ~ 180 MeV (see [O for a review). This smaller value is an artefact of neglecting the radiative 
corrections, which are numerically large. One can hope that these corrections will be reduced significantly 
by taking the ratio of the sum rules. By direct calculations in the limit rrif, — > oo, it has been shown on 
examples of increasing complexity (the Isgur-Wise function WSl E9| and the kinetic energy operator pO|) 
that radiative corrections are indeed very strongly reduced by considering the sum rule ratios, and we 
conform to this practice here. 



pB.- 


■^K* 


pB, 


^<t> 


pB.- 


*K* 



which is not sensitive to the value of the 5-nieson decay constant. Similarly, we have 

0.66 ± 0.09 (45 

and 

pBs^K' 

-^^^ = 0.60 ±0.12. (49) 

The errors given for these ratios are dominated by uncertainties in the 5't/(3)-breaking 
in the wave functions, while the quark-mass dependence is greatly reduced. The small 
value of the last ratio is due to the asymmetry of the K* wave function - the strange 
quark carries a larger fraction of the meson momentum (see Fig. 2). It should be noted 
that errors given here reflect uncertainties in the input parameters, but do not include a 
possible "theoretical" uncertainty of the method itself. We estimate the actual accuracy 
of the calculation to be within 20-30% for the form factors and 10-20% for their ratios. To 
improve the accuracy, it is necessary to calculate corrections to the sum rule (p3|) coming 
from antiquark-gluon-quark components of the vector meson wave functions, S'[/(3)- 
breaking corrections to the wave functions g^ , and the perturbative QCD radiative 
corrections. 

6 The b-quark mass dependence and 
the heavy quark hmit 

In order to facilitate the comparison of the QCD sum rule results with lattice calculations 
0, 1^, we have done the numerical analysis of the sum rules (^) for an arbitrary value of 



the 6-quark mass, assuming the following parametrization suggested by the heavy quark 
effective theory |]52| : 





niB = 


= mb + A, 


A ~ 500-700 MeV 


^o" 


2 


= 2mhUJo , 


cuo^ 1.0-1.2 GeV, 




t = 


= 2mbT , 


r = 0.5-0.8 GeV. 



(50) 

Here A is the binding energy in the limit ruh — > oo, and Uq is the continuum threshold. 
The variable r has the meaning of the non-relativistic Borel parameter |]53|-||57||. 

The results are shown in Fig. 5, where for the reason which is explained below we 
have plotted the value of the form factor multiplied by the meson mass to the power 
3/2. In addition to the curves describing the quark-mass dependence in the heavy-quark 
parametrization (|50[), we also show our result for Ff^^* given in (^61). Furthermore, 
using the mass of the Z^-meson, instead of ms, and the standard values of the parameters 
for D-mesons, s^ = 6 GeV^, rric = 1.4 ± 0.05 GeV, and t = 1.5-3 GeV^ |2|, we obtain 



Ff^^* = 0.85 ±0.13 . (51) 

Although the number given in (^) is not directly physically relevant, it can facilitate the 



comparison with the present lattice data, which are collected mostly in the Z)-meson mass 
range. 



The curves in Fig. 5 disagree slightly with our values for the form factors for B and 



D-mesons. For 5-niesons, the reason is that in (^) we use a rather generous uncertainty 
in the 6-quark mass. For D mesons, the omitted 0{l/m) corrections to the parameters 
given in ( pO]) can already become significant. 

Comparing the values given in ( ^6|) and (|5ll) , we conclude that in the region between 
the charm and beauty masses the form factor scales approximately as Fi ~ l/rriB- How- 
ever, this power behaviour holds only for the sum rule results in the constrained interval of 
quark masses, and should not be confused with the theoretical behaviour in the rrib -^ c)o 
limit, which we are going to discuss now. 



In the heavy 6-quark limit the sum rule in (^3[) involves the integration over a small 
interval of u, of order 1 — m ~ l/fni,. This condition has a simple interpretation: the 
spectator antiquark must recombine with the fast quark coming from the 6-decay. There- 
fore, the form factor is determined by the "tail" of the vector meson wave function cor- 
responding to a strongly asymmetric configuration where almost all the momentum is 



carried by one of the constituents. It has been proven [^, ^ that the behaviour of 
the wave functions in this region at sufficiently large virtualities must coincide with the 
perturbative behaviour, which can be obtained from the conformal properties of the cor- 
responding operators. Quite rigorously, one obtains for the wave functions involved in 
(0): 0x(m) ~ 011 (m) ~ g±^ ~ 0(1 - m), and ^? ~ 0(1) for 1 - m ^ (see (0) and 
(il))- I^ is easy to check that this behaviour is exactly the one needed to ensure that the 
three terms in (E^) are of the same order in the limit rrih — > oo, yielding the form factor 
Fi(0) ~ 0{m^ ). The diagram with the hard gluon exchange produces the same power 



behaviour [0, [T^, and is estimated to be small |T3[. The large factors ~ nih in front 
of the wave functions 0_l and gf in (|23|) are compensated by the smallness of the wave 
functions ~ l/rrih in the relevant region 1 — m ~ l/mf,. 

In order to establish a connection to the heavy-quark expansion of the two-point sum 
rules |53]-[^, we introduce the scale-invariant (up to logarithms) form factor and B- 
decay constant 

Fi(0) = F^{0)■mf, 

Ib = fB-Vm^, (52) 

and rewrite the sum rule ( P^D in terms of the non-relativistic Borel parameter r = t/{2inf,). 
The result reads as: 

/sFi(0)e-^/- = fy j^" due~-'^[-2uj<P'^{l)+mvg'':^{l) - |:mv'(#)'(l)} , (53) 

where the new parameters A and uq are defined in (0). The variable u = mbu/2 has the 
meaning of frequency, and is defined as (p + g)^ — m^ = 2m(,u; (see eq. (p2|)). For the 
wave functions, we obtain, to our accuracy: 

-01(1) = -^0±(1) = 6[l + ai(/x) + (4/5)a2(/x) + (4/3)a3(/x)], 

^?(1) = 3/2, 
-(^f)'(l) = 6, (54) 



where /x ~ 2r [^. In the same hinit, the sum rule for the 5-decay constant becomes 
/le-^/^ = ^ r due--/^u' - {qq), =2r + -^(gagGg)^=2. • (55) 



-m-- 



Using the sum rule in ([55| ) to eliminate the coupling fs from (|53|), we obtain the invariant 
form factor for B -^ K*^ 

Fi(0)~5 (56) 

(with /i^ ~ IGeV^), which translates to the static limit of the decay form factor (at the 
same /x^) 

Fi(0) = Fi{0)/m^J^ ^ 0.4 . (57) 

This number appears to be 30% higher than the result in (46). The scale dependence re- 
minds of possible logarithmic corrections which we have ignored. Note that the standard 
logarithmic factor corresponding to the hybrid anomalous dimension of the interpolat- 
ing current of the i?-meson cancels in the ratio of ( p3D and (0). However, additional 
corrections can appear, similar to the Sudakov-type corrections, which in the context of 
inclusive radiative decays B ^ Xg + •y have been discussed in |TT]] . A related discussion 



of the perturbative corrections produced by the hard gluon exchange can be found in 
||13|| . Since these corrections may exponentiate, the power law in (0) should be taken 
with caution. A study of the perturbative two-loop corrections to the sum rule would be 
welcome, but is beyond the tasks of this paper. 

The contributions of higher twist to the sum rule (p3D do not spoil the existence of the 
heavy-quark expansion. In this case, again, explicit enhancement factors ~ rnl/t ~ 0{inb) 
are compensated by the smallness of the wave functions. However, in contrast with the 
contributions of higher-dimension condensates to the sum rule in (^), the higher- twist 
contributions to (|53| ) generally bring in factors (uj/t), which indicates that the structure 
of the expansion is different. We hope to return to this question in the future. 

Thus, we conclude that the light-cone sum rules, which we suggest in this paper, 
remain well-defined in the heavy 6-quark limit, since all the contributions appear to have 
the same power behaviour at nib -^ oo. The same conclusion has been drawn in the work 
||T9| for the case of semileptonic B — > iiei' decays. 



7 The q^-dependence and the interrelation of 

radiative and semileptonic B-decay form factors 

Up to this point, we have been discussing the radiative decay form factor Fi at g^ = 0. 
Now we are in a position to calculate the g^-dependence. The interest in the g^-dependence 
is due to the fact that within most of the existing approaches the calculations are actually 
done for high values of g^, close to the maximum possible value q'^^^ = (m^ — rriK'Y, 
and then extrapolated to g^ = using plausible assumptions hke the pole-dominance 
approximation. Within the QCD sum rule approach the g^ dependence of form factors 
can be calculated in a wide interval of g^, see e.g. E^, E3. 



The QCD sum rule for the form factor Fi for finite values of q^ is obtained from the 
one given in (^31) by trivial modifications. We obtain 
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r {v)i N , ml + q^ - u^ml + ut („) 1 

The region of applicability of this sum rule is restricted by the requirement that the value 
of q^ — ml be sufficiently less than zero, so that the correlation function (|T^) is evaluated 
in the Euclidean region. In order not to introduce an additional scale in the problem, we 
require that q^ — ml < {p + qY — m^l, which translates to the condition that mf — q^ is 
of order of the typical Borel parameter t ~ 5-8 GeV^. Thus, we end up with an upper 
bound for the applicability of the QCD calculation q^ < 15-17 GeV^, which is a few 
GeV^ below the zero-recoil point. The results are shown in Fig. 6. It is seen that the 
stability of the sum rules becomes worse with the increase of q^. The reason is that the 
omitted contributions of higher twists are suppressed by factors l/{m.l — u{p + g)^) in 
the momentum space, which translates into suppression by powers of ^/{ut) after the 
Borel transformation. At g^ = the average value of u under the integral is of order 0.8, 
but it decreases to ~ 0.6 at q^ c^ 15 GeV^. Thus, in order not to enhance the higher 
twist effects one should somewhat increase the value of the Borel parameter t at high 
g^. This increase, however, potentially comes in conflict with the requirement that the 
Borel parameter t be small enough to allow to separate the i?-meson contribution from 
the continuum. Thus, in general, the light-cone sum rules for large values of g^ are less 
reliable than at g^ = 0. The calculation shown in Fig. 6 was done using the rescaled Borel 
parameter t -^ t/{u){q'^), where (■u)(g^) is the average value of u under the integral. This 
allows for a certain improvement of the stability. 

For the radiative decays discussed here, the sum rules do not contradict the pole-type 
behaviour in the region < g^ < 10 GeV^: 

F^^') = , ^t\ (59) 

1 - gVm^^ig 

with the pole masses of order 3.8-4.2 GeV. A much better quality of the fit is provided, 
however, by the dipole formula 

with mdipoie = 5.6, 5.4, 5.2, and 4.9 GeV for Fig. 6a, b, c, and d, respectively. 

In the case oi B ^ K*'y decays, an important result was derived by Isgur and Wise in 



T^ and by Burdman and Donoghue in ||T3|, where it has been shown that in the heavy 



6-quark limit an exact relation exists between the form factor Fi and the semileptonic 



form factors defined by: 
{K*{p',X)\-s^,{l-^,)b\B{p)) 



2V{q^ 



-fiupa^ 



*(A).pPp'-_,(^^+^^,)A^(g2)g*(A)^_ 



where the elhpses denote terms proportional to [p + p')^ or q^. According to |]12|, |13|: 

Viq' 



(61) 
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2m 
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(62) 



B 



as far as the value of g^ is sufficiently close to the point ql 



{ms - mK*f 
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GeV^. Similar relations can of course be derived for other radiative decays. The central 
problem in using predictions based on the heavy quark effective theory to the rare B- 
decays is whether this relation is satisfied for small values of g^, as suggested in |[T^, [T^, 
and what is the size of the power l/rrib corrections. Our main result in this section is that 
the relation in (|62|) is fulfilled in the QCD sum rule approach discussed here to a good 
accuracy for finite 6-quark masses and for all momentum transfers. 

The sum rules for the semileptonic form factors entering (|62D can be derived in a 
similar way, repeating the steps described in section 3. A simple calculation yields 
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(64) 



Note that the three sum rules in ([58|), (| 
specified in section 3. 

It is now easy to check the Isgur-Wise relation by forming the ratio 

i?(g2)=Fi(g2)/Fr(g'), 



and (I6l) involve the same wave functions as 



(65) 



where -Fi(g^) is the form factor of the radiative decay according to (|58|) , and FY^{q^) is 
the particular combination of the semileptonic form factors, which enters on the r.h.s. of 

) and (|6^ . The ratio (|65D is 



(p^), and which is evaluated by using to the sum rules ( 
shown as a function of g^ in Fig. 7 for the two decays B -^ K*'~f and B -^ p'j, and is 
practically the same for the two remaining decays Bg — > 07 and Bs — ^ K*'~f. We find 
that the relation between radiative and semileptonic form factors is satisfied with a high 



accuracy for all momentum transfers. For the physically relevant point g^ = a typical 
value of the ratio R is about 0.92-0.95. These predictions should be rather reliable, since 
many of the uncertainties in the individual decay rates are greatly reduced in the ratios. 
Since the semileptonic decay B — > K* + i + Ui does not occur in nature, tests of the 
Isgur-Wise relation must take into account the correction for the SU{3) breaking between 
B -^ K* and B —>■ p semileptonic form factors. From our sum rule calculation we get 



\B^K* 



A(o)^^7Ai(o)^ 
\/(o)^^7r(o)^^^* 



0.76 ±0.05, 
0.73 ±0.05. 



(66) 



As a by-product of our analysis, we have derived the q^ dependence of the semileptonic 
decay form factors, V and Ai, which is shown in Fig. 8 for the transition B ^ p + C. + v^. 
In Table 1 the values for V{q^ = 0) and Ailq"^ = 0) are compared with the corresponding 
results from other approaches @^[§31- ^^ ^^^ ^ satisfactory agreement with the quark 
models, while our results are somewhat lower than previous QCD sum rule estimates, 
which were done using the traditional approach. In agreement with [^ we obtain a 
steeper increase of the form factor Vlq"^) with g^, compared to y4i(g^), although values of 
the slopes come out to be different. 

Predictions of the traditional sum rules should be more reliable at the values ml — q^ = 
0{mb) than at g^ = 0, since in this region they do not suffer from large contributions 
of the operators of high dimension, see the detailed discussion in the appendix. In the 
light-cone sum rules the situation is opposite, so that our results are to a large extent 
complementary to the calculation in []58| , |60[] . It is worth while to note that the difference 
in predictions is minimal at intermediate q^ ~ 5-10 GeV^, i.e. exactly in the region where 
both the approaches should work well. 



Table 1 : Form factors for the semileptonic decay B ^ p + i + i>£ 



Reference 


Ai(0)«^^ 


y(o)^^^ 


This paper 


0.24 ±0.04 


0.28 ±0.06 


58" 


0.5 ±0.1 


0.6 ±0.2 


59'^ 


0.35 ±0.16 


0.47 ±0.14 


3C ' 


0.28 


0.33 


3^5 


0.05 


0.27 


0^ 


0.21 


1.04 



a: QCD sum rules 

b: Quark models 

c: HQET + chiral perturbation theory 



8 Summary and Conclusions 

Motivated by the discovery of the electromagnetic penguins through the decay B -^ 
K* +7, we have presented an alternative method of calculating the transition form factors 
in this and related decays using the approach of QCD sum rules on the light cone. Our 
calculations give F^-*^* = 0.32 ± 0.05; combined with the estimates of the inclusive 
branching ratio BR{B ^ X^ + 7) = (3.0 ± 1.2) x 10"^ in the Standard Model [|], this 
yields BR{B —>■ K* +7) = (4.8 ± 1.5) x 10~^, which is in reasonably good agreement with 
the observed branching ratio (4.5 ±1.5 ±0.9) x 10~^ measured by the CLEO collaboration. 

We use the same method to calculate the related form factors in the decays B^^a -^ 
p + •y, Bd ^ uj + -y, Bs ^ (f) + •y and Bg -^ K* + 7. These decays can be related to the 
already discussed and measured decay Bu^ — ^ K* + 7. To this end, we have evaluated 
in detail the numerical effect of the S'f/(3)-breaking in the form factors. This fixes the 
decay width T{Bs — ^ + 7), which involves the same CKM matrix element \Vts\ as the 
decay width r{Bu^d -^ K* + 7). Since the form factors for the two decay modes are 
comparable in size, see Fig. 3 and eq. (|46|), and since the lifetimes of the B^iBd and Bg 
mesons are expected to be very similar — an expectation which is consistent with present 
experiments — we estimate that BR{Bs — i> + 7) ~ BR{Bu^d -^ K* + 7). 

More interesting, however, are the CKM-suppressed decay modes discussed above. 
A measurement of both the CKM-allowed and CKM-suppressed radiative decays B -^ 
V + y would give direct information on the CKM matrix element ratio |l^d|/|Vis|. The 
form factors needed for such determination are plotted in Fig. 4, showing the ratios 
of Fj^^^/Ff ^^*, Ff'^^'/Ff'^'^ and pBs^K* 1 pB^A^K* ^ ^^^ numerical results for the 
form factor ratios are given in eqs. (^7|) -(^9l). The ratios of the branching ratios BR{B — >■ 
Lu + y)/BR{B -^ K* +y) and BR{Bs ^ K* + y)/BR{B -^ K* + y) as functions 
of |V^d|/|V^s| are shown in Fig. 9 {B = B^ or -B^.). As already stressed, these ratios 
are practically independent of rrit and of the QCD scale parameter fi. The dependence 
of the ratios on the various parameters endemic to the QCD sum rule approach also 
becomes rather mild, as opposed to the predictions for the partial decay widths themselves. 
Therefore, the two curves shown in Fig. 9 should give a fair representation of the residual 
theoretical error, and they allow a rather clean determination of the ratio of the CKM 
matrix elements |Vtd|/|Vfs| as and when the required data for the ratios of branching ratios 
become available. 

It should be pointed out here that similar relations have been derived for the ratios 
of the inclusive branching ratios BR{B -^ X^ + y)/ BR{B ^ X^ + 7) in [Q, and for the 
ratios of the B^-B^ mixing parameters Xd and Xs |0. In particular, the Standard Model 
predicts: 

Xd _ ^b,Mb, {Pb.Bb, 



Xs tbMb. {PbBb. 



Vu'' 



V, 



ts 



(67) 



where Bb^ and Bb^ are the hadronic bag constants. Since all dependence on the t-quark 
mass and the QCD renormalization of the box amplitudes drops out, we are left with 
the square of the ratio of CKM matrix elements, multiplied by a factor that reflects 
S'[/(3)flavour-breaking effects. With the known values of Xd and BR{B — > K* + 7), a 
measurement of Xg and of the CKM-suppressed radiative decays would provide valuable 



quantitative tests on the CKM matrix and, in particular, its unitarity. 

In order to obtain an estimate of the absolute branching ratios for the CKM-suppressed 
decays, we use the constraints on the matrix element \Vtd\ that already exist from the 
measurement of the mixing ratio Xd- The present value Xd = 0.71 ± 0.07 [^ yields the 
model dependent bound 0.005 < \Vtd\ < 0.012 |^. This in turn gives 



0.10 < \-f\ < 0.33, (6^ 

where we have used the CKM unitarity constraint and recent data giving |Vts| — |Vcfe| - 



0.042 ± 0.005 |33|. We can combine this bound with Fig. 9 to predict: 



BR(B ^ cj + 7) 

with the ratio BR{B -^ p + j)/BR{B — > K* + 7) bounded likewise. While this gives 
a ball-park estimate for the branching ratios in question, it is clear from this numerical 
exercise that a measurement of the above ratio within a factor 2 would enormously reduce 
the present (model-dependent) uncertainty on the CKM matrix element ratio |V^(i|/|V^s|- 
The main theoretical content of this paper is in pointing out and providing general 
arguments for the light-cone QCD sum rule approach to be a more adequate tool for 
the study of heavy meson decay form factors compared to the traditional approach. The 
principal input in the light-cone sum rules are the wave functions of vector mesons, and 
the accuracy of our calculation can be substantially improved by taking into account the 
wave functions of twist 3 and twist 4. For the case of the 7r-meson, a complete set of wave 
functions to this accuracy is given in ^^. For vector mesons, this problem is not solved 



yet, and is of acute practical interest. Some new results concerning the interrelation of 
the transverse and longitudinal wave functions of vector mesons are given in section 4. 

We have also studied the heavy-quark mass dependence of Fi in the region rric < 
rriQ < rrih. This is useful for comparison with lattice-QCD results where the heavy quark 
systems at present are calculated in the vicinity of rric and extrapolated to mf,. We find 
that the transition form factor Fi scales approximately in this region as Fi ~ l/rrib, and 
has a value Fi = 0.85 ± 0.13 for rriB equal to the D-meson mass, i.e. for the (unphysical) 
process D -^ K* + 7. 

Apart from the phenomenologically interesting value of the form factor F^~*^ at the 
on-shell point g^ = 0, we have investigated the behaviour of Fi as a function of g^. The 
case g^ 7^ contributes to the decays B ^ V + £~^i~ via off-shell photons. In addition, 
we have derived the light-cone QCD sum rules for the form factors ^1(5'^) and V^(q'^), 
describing the semileptonic decay B ^ p + iui. It will be interesting to compare them 
with data on the semileptonic decay B ^ p + iu^, as and when these become available. 
Another use of this sum rules is to check the relation between the radiative decay form 
factor Ff^^ and the semileptonic form factors in B —>■ p + iui, which was derived 
in |12, 13] using the heavy-quark symmetry and for g^ near the kinematic point g^ax- 
We have shown that the ratio Fi(g^)/F/^(g^) evaluated in the QCD sum rule approach 
remains very close to unity in the complete kinematic range of g^ which is available in the 
CKM-suppressed semileptonic decays. 



In summary, we have given a number of predictions here in the crucial area of ex- 
clusive radiative and semileptonic 5-decays, involving the CKM-suppressed and CKM- 
allowed transitions using QCD sum rules on the light-cone. These predictions are specific 
to this approach and it will be instructive to confront them with data. The existing 
data on BR{B -^ K* + 7) are reasonably well explained, and we argue that our results 
allow to quantify the ratios of CKM-allowed and CKM-suppressed exclusive radiative 
5-decays. This would prove useful in determining the CKM ratio |l^td|/|^s| in a fairly 
model-independent way. 

Acknowledgements 
We acknowledge helpful discussions with David Cassel, Michael Danilov, Thomas Mannel, 
Stephan Narison and Matthias Neubert. One of us (V.M.B.) would like to thank John 
Ellis and the CERN TH division for their hospitality. 



Appendix 



Three-point sum rule 

In this appendix we compare our approach, for the particular case oi B ^ K*'y decays, 
with the results of the standard QCD sum rule analysis |§|, ||, 0]- Here, the starting point 
is the operator product expansion for the three-point function involving the penguin 
operator and the interpolating currents for both the i?-meson and the K*-TLneson 

z^ fdxdye'P^+'^y{0\T{[dj^s]ix)[sa^,qMy)[bil5dmm (A.l) 



in which the contribution of interest is identified with the one having poles at p^ = m\* 
and (p + qY = fn^B- The standard treatment, details of which we are not going to present 
here, leads to the following sum rule: 

rrib ^ 

/ ° dsi e-^i/*i / ° ds2e-^'^-<^''^e{s2 - si - m^)p(si, Ss) - (m, + mh){dd) 

Jo Jmi 



47r^ Jo «,„.j 



^-i,a,G^ °\ ' " H- "j; ■' ,^-^.-^, (A.2) 



where p{si, S2) is the spectral density of the triangle diagram in Fig. 10a 

/ N ^bSi , S2-Si-ml 
p{si,S2) = - -r + mfcm,— -^, (A.3) 

[S2 — Si) [S2 — Si) 

and we have calculated also the contribution of the quark condensate in Fig. 10b, and of 
the mixed condensate in Fig. 10c (a typical graph is shown). The contributions of the 
gluon condensate and of the four-quark condensate are negligible [0]. In contrast with 
ref. p|, we give the answer with ©(m^) corrections included. To the accuracy that these 
corrections are put to zero, our sum rule agrees with the one given in [^, but for the sign 
in front of the third term in the contribution of the mixed condensate. Numerically, the 
difference is negligible. Note that in contrast to the light-cone sum rules we have now two 
Borel parameters, ti for the i^*-meson and ^2 for the 5-meson, which should be taken 
two times larger than the Borel parameters in the corresponding two-point sum rules (cf. 
p5|). The quantity Sq = 1.7 GeV^ is the continuum threshold in the sum rules for the 
i^*-meson, and rUs is the mass of the strange quark, which is of order rus = 150 MeV. 
The other entries have been specified in the text. 

A numerical treatment of the sum rule in (|A.2| ) yields values of the form factor of 
order 

Ff^^*(0) -0.5-0.6 (A.4) 



which is in agreement with the analysis in @, 0. The lower value ~ 0.38 quoted in [0] was 
obtained by the rescaling of a similar value ~ 0.55 from the scale of the typical Borel pa- 
rameter to the scale fi = rrih, using the large anomalous dimension of the penguin operator 
(I). We have not understood the reasoning for such a rescaling, since the renormalization 
group treatment of the effective Hamiltonian in @ makes sense at scales /i > mh only. 
A somewhat smaller value obtained in ref. 0] is due to a smaller value of the quark 
condensate and due to a higher value fs = 180 MeV used in p. 

The value in ( [A.4|) appears to be considerably larger than our result in (|46|). We are 



going to claim, however, that the sum rule in ( [A.2|) is ill-behaved in the limit of heavy 
6-quark, and is less reliable. 

An inspection of the sum rule in ( A. 2 ) shows that various terms in it have a different 



behaviour in the limit rrih -^ oo. One easily finds that the perturbation theory contribution 
to ( |A.2| ) yields the form factor -Fi(O) ~ 0(m^ ), the quark condensate produces a 
contribution of order m^ , and the mixed condensate contribution is of order m^ . Thus, 



the sum rule ( [A.2| ) blows up in the limit mb -^ oo. 

This problem has actually been known for a long time, and was originally found in 



the QCD sum rule calculations of the pion form factor [^, |2^. In technical terms, 
the problem appears because the operator product expansion for three-point functions 
involves the dimensionless parameter (??i^ — q'^)/t, increasing powers of which multiply 
the contributions of local operators of higher dimensions. In the case of light quark 
systems, m^ = and the Borel parameter is of order 1 GeV^. Thus the expansion breaks 
down at sufficiently large values of q^. In the case of heavy-light mesons t ~ Olrrih), and 
for g^ = the operator product expansion involves increasing powers of the heavy-quark 
mass. In the present case, we observe a considerable enhancement of the contribution of 
the mixed condensate, and conclude that higher-order corrections, e.g. proportional to 
the dimension-7 condensate (ipg'^G'^il)), can influence the result significantly. It has been 
claimed recently ||6^ that this condensate is much larger than its factorized value. 



In physical terms, the reason of the difficulty with the traditional three-point sum rules 
in the rrih — >■ oo limit is that the operator product expansion does not take into account 
the effect of the finite correlation length between the quarks in the physical vacuum. 



A possible remedy, suggested in [^, is to use the expansion in non-local condensates, 
which take into account this effect in a mo del- dependent way, see e.g. [^. Following 
Radyushkin, we write down the non-local quark condensate in the Borel representation: 

{d{x)diO)) = {dd) dve'''' l^f{v). (A.5) 



Moments of the function /(z/) are determined by vacuum expectation values of local 
operators 



rfz//(z/) 


oo ][ _ 



/•oo I _ _ 

/ dvvf{v) = -{dagGd)/{dd) c^ 0.2 GeV\ (A.6) 

Jq 4 

and similar relations for higher moments. An effect in the sum rule ( |A.2D is that in the 



term proportional to the quark condensate one should make a replacement 

Idd) -^ {dd) rrfz//(z/)exp|-^^^| (A.7) 

Jo [ h t — u j 

where t = ^1^2/(^1+^2), and discard the contribution of the mixed quark-gluon condensate, 
since its contribution is partly included in (|A.5| ). In the limit of the heavy 6-quark mass 
the integral over u in ( [A.7| ) is dominated by contributions of small z/; thus the behaviour of 
/(z/) at z/ — > is of crucial importance. Note that no information about this behaviour can 
be extracted from known values of a few first moments. In a confining theory such as QCD 
one should expect that correlation functions decrease exponentially at large distances in 
Euclidean space, 

{d{x)d{0)) "' ~"°° exp{-a^^^} , (A.8) 

where a is the correlation length. Existence of a finite correlation length in the physical 
vacuum is a part of a common wisdom about confinement and a starting point of many 
existing models of the vacuum (see e.g. [Q). It is easy to show that the finite correlation 
length implies quite generally the exponential behaviour of the function /(z/) at small z/: 

/(z/) ~ e^"'/^ (A.9) 

Inserting this behaviour in (|A.7| ) and taking the limit nib ^ 00 we obtain to the expo- 



nential accuracy the contribution of the quark condensate to be 



~ exp{-2amfc/VM2} . (A.IO) 

Since ti = 0(1) and ^2 = 0{mb), this result means that the contribution of the quark 
condensate decreases as exp{ — y^m^}- Hence, the perturbative contribution dominates 
the sum rule in this limit, yielding the same behaviour -Fi(O) ~ m^ which is intrinsic to 
the light-cone sum rules in the text. Note that this result is a quite general consequence 
of the existence of a finite correlation length in the QCD vacuum. It is actually not 
surprising, since the power behaviour of exclusive process with large momentum transfers 



2l|] is given by perturbative diagrams. In the present case, the diagrams with hard 



gluon exchange yield the same power behaviour as the ones considered here, see |19 |. 

To summarize, we argue that the light-cone sum rule written in (^3D remains well- 
defined in the limit of the heavy 6-quark mass, similar to the traditional QCD sum rules 
for the two-point functions, which can be written directly in the heavy-quark effective 
theory |^|. The heavy-quark limit does not exist, however, in the standard sum rules 



for three-point functions, as the one in ( [A.2|) . This property of the light-cone sum rules 
argues in favour of our approach. 
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Figure captions 

Fig. 1 The leading contribution (a) and the gluon correction (b) to the correlation func- 
tion in (|T5|). 

Fig. 2 The leading-twist wave function 0_l (p!9D from the sum rule analysis in [^ at the 
scale /i^ = 5 GeV^ for (a) p- (solid line) and 0- (dashed dotted hne) mesons and (b) 
for the ii'*-meson. The wave function for K* is also shown at the low scale /x^ = 1 
GeV^ (dashed line). The dotted lines correspond to the asymptotic wave function 
01" = 6m(1-m). 



Fig. 3 Stability plots for the sum rule in (p3D as a function of the Borel parameter for 



B — > -ft'*7, B — > p7, Bs -^ 07, and Bs — > -ft'*7, see a, b, c and d, respectively. 
Dashed lines correspond to Sq = 35 GeV^ and dotted hues to Sq = 33 GeV^ (with 
Sq" given in (|4^)). For each value of Sq curves are shown for nif, = 4.6 GeV (upper 
curve), rrih = 4.7 GeV (middle) and rrih = 4.8 GeV (lowest curve). 

Fig. 4 Ratios of the form factors i^i(O) for the processes (a) -B — *> p7 and B -^ K*'y, 

(b) Bs -^ K*^ and Bg -^ 07, and (c) Bg -^ K*'y and B -^ K*'y. Dashed lines 
correspond to Sq = 35 GeV^ and dotted lines to s^ = 33 GeV^ (with Sq" given in 
(^2l)). The corresponding upper curves are for m^ = 4.6 GeV and the lower ones 
for nib = 4.8 GeV. 

Fig. 5 Dependence of the form factor -Fi(O) on the quark mass. The various parameters 
are scaled according to (^D|). Solid lines correspond to r = 1.0 GeV and the dotted 
lines to r = 0.6 GeV. The corresponding upper curves are for A = 600 MeV with 
uo = 1.2 GeV, and the lower ones for A = 500 MeV with ujq = 1.0 GeV. For 
comparison the results corresponding to the parameters discussed in section 5 are 
indicated by vertical bars. 

Fig. 6 Momentum dependence of the form factors Fi{q'^) for (a) B —>■ K*^, (b) B — »• pj, 

(c) Bg — >■ 07, and (d) Bg — > K*'j. Dotted lines correspond to t = 5GeV^/(u) 
and dashed lines to t = 8 GeV'^ /{u). The corresponding upper curves are for rrib = 
4.6 GeV and the lower ones for mi, = 4.8 GeV. 

Fig. 7 Momentum dependence of the ratio of radiative and semileptonic form factors (a) 
B -^ K*'j and (b) B -^ p'j. Dotted lines correspond to t = 5 GeV'^/{u) and dashed 
lines to t = 8 GeV^/(n). The corresponding upper curves are for nib = 4.8 GeV 
and the lower ones for nib = 4.6 GeV. 

Fig. 8 Momentum dependence of the form factors ^1(5^) (a) and Vlq"^) (b) for B -^ p'y. 
Dotted lines correspond to t = 5 GeV'^ / (u) and dashed lines to t = 8 GeV'^/{u). 
The corresponding upper curves are for nib = 4.6 GeV and the lower ones for 
mb = 4.8 GeV. 

Fig. 9 Ratio of the branching ratios BR{B -^ uj'~f)/BR{B —>■ K*'~f), where B = Bu 01 
Bd, and BR{Bg -^ K*j)/BR{B -^ K*j), as a function of the CKM matrix elements 
|^d|/|Vfs|. The two lines correspond to the errors given in (^7]). 



Fig. 10 Contributions of perturbation theory and of vacuum condensates to the sum 
rule in (Q. 

Fig. 11 Y Stabihty plot for the sum rule in (|A.2D 
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